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(57) ABSTRACT 

Computational aeroelastic analyses typically use a math- 
ematical model for the structural modes of a flexible structure 
and a nonlinear aerodynamic model that can generate a plu- 
rality of unsteady aerodynamic responses based on the struc- 
tural modes for conditions defining an aerodynamic condition 
of the flexible structure. In the present invention, a linear 
state-space model is generated using a single execution of the 
nonlinear aerodynamic model for all of the structural modes 
where a family of orthogonal functions is used as the inputs. 
Then, static and dynamic aeroelastic solutions are generated 
using computational interaction between the mathematical 
model and the linear state-space model for a plurality of 
periodic points in time. 

27 Claims, 6 Drawing Sheets 
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METHOD OF PERFORMING 
COMPUTATIONAL AEROELASTIC 
ANALYSES 

Pursuant to 35 U.S.C. §119, the benefit of priority from 5 
provisional application 60/908,769, with a filing date of Mar. 

29, 2007, is claimed for this non-provisional application. 

ORIGIN OF THE INVENTION 

to 

This invention was made by an employee of the United 
States Government and may be manufactured and used by or 
for the Government of the United States of America for gov- 
ernmental purposes without the payment of any royalties 
thereon or therefor. 15 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to computational aeroelastic analy- 20 
ses. More specifically, the invention is a method for perform- 
ing computational aeroelastic analyses using linear state- 
space models. 

2. Description of the Related Art 

When designing aerodynamic structures, it is important to 25 
understand if and under what aerodynamic conditions the 
structure will be stable and unstable. This is especially true if 
a structure is inherently flexible (e.g., a wing) such that 
unstable operation thereof leads to oscillations of the struc- 
ture until failure occurs. Accordingly, the field of aeroelastics 30 
examines the interaction between a flexible structure and the 
unsteady aerodynamics generated by the structure. Since 
structure design generally includes performance predictions, 
it is necessary to predict the aeroelasticity of a particular 
design. To do this, computational aeroelastic methods are 35 
used to numerically simulate an aeroelastic process using 
computational techniques. In general, the computational 
techniques include the use of a mathematical model of the 
structure and an aerodynamic model of aerodynamic condi- 
tions of the structure. 40 

Referring now to FIG. 1, traditional computational 
aeroelastic analyses use Computational Fluid Dynamics 
(CFD) codes that require the coupled feedback interaction of 
a linear structural model 10 and a nonlinear aerodynamic 
model 20 (i.e., the CFD code). This coupled feedback inter- 45 
action consists of information 12 (e.g., structural modes such 
as displacements or deflections in the case of flexible struc- 
tures) from structural model 10 being passed to nonlinear 
aerodynamic model 20. For a flexible aerodynamic structure 
such as a wing, information 12 from structural model 10 50 
consists of the physical displacements of the structure due to 
an aerodynamic force. This physical displacement informa- 
tion is supplied to nonlinear aerodynamic model 20 in order to 
compute an aerodynamic response or force 22 induced by this 
structural displacement. The computed aerodynamic 55 
response 22 is then fed back to structural model 10 to compute 
new displacement information 12. This exchange of informa- 
tion repeats at each time step of a numerical solution and this 
coupled feedback interaction takes place wholly within the 
aeroelastic CFD code. 60 

Nonlinear aerodynamic model 20 is significantly more 
complicated than linear structural model 10. Therefore, the 
computational efficiency of the coupled feedback interaction 
process is driven by the large computational expense of the 
nonlinear aerodynamic system defined by the CFD code. As 65 
a result, traditional computational aeroelastic analyses are 
very time-consuming and computationally expensive. The 
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output from the repetitive feedback process described above 
are time histories of the response of the structure at a given 
flight condition. The output time histories indicate the level of 
stability (or instability) of the structure’s configuration being 
analyzed at a given flight condition. Stability (instability) is 
determined by the convergence (divergence) of the time his- 
tories. By changing the flight condition and repeating this 
analysis as necessary, the region of the flight envelope where 
the configuration is safe (stable) or unsafe (unstable) is 
defined. Clearly, computation of the stability/instability 
boundary (i.e., also known as flutter) becomes very expensive 
and time-consuming since the above-mentioned process must 
be repeated at several conditions. 

In addition, the output time histories are in a form that is not 
suitable for use by other structure-design disciplines such as 
controls or optimization. These disciplines have a very spe- 
cific requirement for the types of mathematical models that 
can be analyzed. However, the output time hi stories from the 
traditional computational aeroelastic analyses are not in that 
form. Thus, the value of the information generated by tradi- 
tional computational aeroelastic analyses is limited in that it 
cannot be readily utilized by other disciplines involved in the 
overall vehicle design process. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the present invention to 
provide a method of performing computational aeroelastic 
analyses. 

Another object of the present invention is to provide a 
computationally-efficient method of performing aeroelastic 
analyses that reduces the amount of time it takes to generate 
an aeroelastic solution. 

Still another object of the present invention is to provide a 
method of performing aeroelastic analyses in order to gener- 
ate an aeroelastic solution that is in a form suitable for use by 
other disciplines in the aerodynamic structure design process. 

Other objects and advantages of the present invention will 
become more obvious hereinafter in the specification and 
drawings. 

In accordance with the present invention, a method is pro- 
vided for performing aeroelastic analyses. The method uses a 
mathematical model for a plurality of structural modes of a 
flexible structure, and a nonlinear aerodynamic model that 
can generate a plurality of unsteady aerodynamic responses 
based on the structural modes for conditions defining an 
aerodynamic condition of the flexible structure. A linear 
state-space model is generated using a single execution of the 
nonlinear aerodynamic model for all of the structural modes. 
This single execution of the nonlinear aerodynamic model is 
carried out using a family of orthogonal functions as inputs 
thereto. Then, static and dynamic aeroelastic solutions are 
generated using computational interaction between the math- 
ematical model and the linear state-space model for a plural- 
ity of periodic points in time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a simplified graphic depiction of a conventional 
computational aeroelastic process for a flexible structure; 

FIG. 2 is a graphic depiction of a method of performing 
computational aeroelastic analyses for a flexible structure in 
accordance with an embodiment of the present invention; 

FIG. 3 illustrates a set of orthogonal block pulse functions; 

FIG. 4 illustrates a set of orthogonal Haar functions; 

FIG. 5 illustrates a set of orthogonal Walsh functions; and 
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FIG. 6 is a graphic depiction of a method of performing 
computational aeroelastic analyses in accordance with 
another embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Referring again to the drawings and more particularly to 
FIG. 2, a method of performing computational aeroelastic 
analyses in accordance with an embodiment of the present 
invention is illustrated graphically. The attributes of this novel 
method that are the same as the traditional method (FIG. 1) 
will be referred to using common reference numerals. This 
new method can be applied to any flexible structure that can 
be modeled in terms of its various structural modes or dis- 
placement characteristics. For example, the present invention 
can be used to perform computational aeroelastic analyses on 
a variety of flexible aerodynamic structures such as wings, 
fuselages, and tail components. 

Similar to traditional computational aeroelastic analyses 
methods, the present invention uses a structural model 1 0 that 
is a mathematical model describing the various structural 
modes for a flexible structure of interest. More specifically, 
structural model 10 provides an aeroelastic model that can be 
used within and by the analyses’ aerodynamic model. For 
example, if the aerodynamic model is a “computational fluid 
dynamics” (CFD) code, structural model 10 will include a 
CFD grid and structural mode shape information as is known 
in the art. During any given time step of the process illustrated 
in FIG. 2, structural model 10 presents displacement infor- 
mation 12 that describes a structural mode shape in terms of 
the CFD grid. 

In accordance with the present invention, all of the struc- 
tural modes defined by structural model 10 are initially pro- 
vided to a nonlinear aerodynamic model 20 that can be (and 
typically is) a CFD code. A well known and proven example 
of a CFD code is the CFL3D code developed by the National 
Aeronautics and Space Administration (NASA). This CFD 
code is available from NASA Langley Research Center, 
Hampton, Va. 

The present invention requires and performs only a single 
execution of model 20 for all of the structural modes of 
structural model 10. All of the structural modes are excited 
simultaneously using inputs that are defined by a family of 
orthogonal functions 30. From a mathematical perspective, 
two functions are considered to be orthogonal when the wave- 
forms or signals associated therewith are as different or inde- 
pendent from one another as possible. Three types of func- 
tions that can be used in the present invention are the known 
block pulse functions, Haar functions and Walsh functions. 
By way of example, FIG. 3 illustrates a set of four block pulse 
functions (that would be used to solve a four input problem), 
FIG. 4 illustrates a set of four Haar functions, and FIG. 5 
illustrates a set of four Walsh functions. 

The single execution of model 20 for all of the structural 
modes of model 10 using orthogonal inputs 30 produces 
unsteady aerodynamic responses 32 for all of the structural 
modes simultaneously. Responses 32 are provided to a con- 
version process 40 that transforms responses 32 to a linear 
model 50 in state-space form. The state-space form is a well- 
known reduced-order model form that is also used by other 
aerodynamic design disciplines to include the controls and 
optimization disciplines. In the present invention, linear state- 
space model 50 is a simplified mathematical model that rep- 
resents unsteady aerodynamic responses 32 extracted from 
nonlinear aerodynamic model 20. In other words, state-space 
model 50 has the characteristics of nonlinear aerodynamic 
model 20. However, by its very nature, linear state-space 
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model 50 generated by conversion process 40 is significantly 
less complex than nonlinear aerodynamic model 20. 

Conversion process 40 can be implemented using a variety 
of routines known to those skilled in the art. For example 
5 conversion process 40 can be realized using the Eigensystem 
Realization Algorithm (ERA). See Juang et al., “An Eigen- 
system Realization Algorithm for Model Parameter Identifi- 
cation and Model Reduction,” Journal of Guidance, Control, 
and Dynamics, Vol. 8, 1985, pp. 620-627. In terms of execut- 
10 able computer code, the ERA is a subroutine of the M ATL AB 
software package available commercially from The Math 
Works, Inc., Natick, Mass. Briefly, the ERA is a routine used 
to create a linear math model of a system if the inputs and 
corresponding outputs of the system are known. 

15 The present invention next applies the traditional compu- 
tational feedback interaction process using structural model 
10 and linear state-space model 50 in order to generate static 
and dynamic aeroelastic solutions. A well known platform/ 
environment for carrying out this process is the MATLAB/ 
20 SIMULINK environment, also available commercially from 
The Math Works, Inc. It is to be understood that other soft- 
ware tools could also be used to carry out the feedback pro- 
cess without departing from the scope of the present inven- 
tion. Since the MATLAB/SIMULINK environment uses the 
25 mathematical state-space definition of systems to perform 
simulations and analyses, it is a good choice to use with linear 
state-space model 50. This feedback process is illustrated in 
FIG. 2 by the outermost “circle” of events. As would be 
understood in the art, this feedback process runs over a plu- 
30 rality of time steps to arrive at aeroelastic solutions. 

Referring still to FIG. 2, structural model 10 should be in a 
state-space form in order to be compatible with linear state- 
space model 50. If structural model 10 is not in state-space 
form, it is necessary to apply a state-space form converter 11 
35 to structural model 10. Methods for putting structural model 
10 in a state-space form are well known in the art. The result- 
ing displacement information 12' in state-space form is 
passed to linear state-space model 50. 

The linear nature of both displacement information 121 
40 and model 50 allows the resulting aeroelastic responses 52 to 
be generated much faster than an execution of nonlinear aero- 
dynamic model 20. Aeroelastic responses 52 are then passed 
to structural model 10 where the process repeats itself. As 
mentioned above, the process repeats for a number of time 
45 steps until the resulting static and dynamic aeroelastic solu- 
tions indicate system stability or instability. Solution times 
are greatly reduced in the present invention as substantial 
computational efficiencies result from using linear state- 
space model 50. At the same time, since linear state-space 
50 model 50 captures the essence of the computationally-inefli- 
cient nonlinear aerodynamic model 20, results generated by 
the present invention are similar to those generated by the 
traditional computational feedback interaction process that 
relies solely on nonlinear aerodynamic model 20. Further- 
55 more, the resulting aeroelastic solutions are in state-space 
form and are, therefore, in a form that is readily usable by 
other aerodynamic control disciplines. 

The computational speed/efficiency of the present inven- 
tion can be further enhanced by the provision of predeter- 
60 mined aerodynamic offsets 54. Offsets 54 are predetermined 
values indicative of an aerodynamic load that would result if 
the structure of interest were a rigid body. In other words, 
offsets 54 essentially define the aerodynamic load on a rigid 
or non-flexible embodiment of the structure. The determina- 
65 tion of such offsets/values is well understood in the art. By 
starting the feedback process with rigid aerodynamic load 
values, the actual aeroelastic (i.e., flexible structure) static 
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response solution will be representative of the true physical 
deformation of the structure. Often, in an aeroelastic analysis 
process, either static aeroelastic solutions or static plus 
dynamic aeroelastic solutions are desired. Due to the innova- 
tion of using the rigid aerodynamic offsets, the present inven- 
tion can be used to provide either static aeroelastic solutions 
or static plus dynamic aeroelastic solutions. This is a new 
capability for aeroelastic analyses using state-space models. 

The present invention also provides the basis for additional 
computational advantages. In particular, the present invention 
can be expanded to provide what are known as matched-point 
solutions, i.e., new aeroelastic solutions for changes in veloc- 
ity. To do this, the above-described process is repeated for 
different time step spacings or periods. More specifically and 
with reference to FIG. 6, after the above-described process is 
used to generate aeroelastic solutions for a given time step 
period At l3 the time step period At is changed at 60 and a new 
aeroelastic solution set is generated by the outermost “circle” 
of events described above. Note that the time step must be 
changed for both linear state-space model 50 and the state- 
space form of structural model 10 since both models must run 
in the same time step. Thus, the present invention can also be 
used to generate matched-point solution sets with only a 
single execution of nonlinear aerodynamic model 20. 

Additional discussion of the present invention is provided 
in Walter A. Silva, “Simultaneous Excitation of Multiple- 
Input Multiple-Output CFD-Based Unsteady Aerodynamic 
Systems,” 48* AIAA/ASME/ASCE/AHS/ASC Structures, 
Structural Dynamics, and Materials Conference, AIAA 2007- 
1988; and Walter A. Silva, “Recent Enhancements to the 
Development of CFD-Based Aeroelastic Reduced-Order 
Models,” 48'* AIAA/ASME/ASCE/AHS/ASC Structures, 
Structural Dynamics, and materials Conference, AIAA 2007- 
205 1 , both herein incorporated by reference in their entirety. 

The advantages of the present invention are numerous. The 
time required to generate static and dynamic aeroelastic solu- 
tions to computational aeroelastic analyses is drastically 
reduced from hours/days to minutes/hours. The generated 
aeroelastic solutions are in a state-space form that is compat- 
ible with the form of data used by other aerodynamic design 
disciplines, thereby making the overall aerodynamic design 
process more efficient when the present invention is used 
therewith. 

Although the invention has been described relative to a 
specific embodiment thereof, there are numerous variations 
and modifications that will be readily apparent to those 
skilled in the art in light of the above teachings. It is therefore 
to be understood that, within the scope of the appended 
claims, the invention may be practiced other than as specifi- 
cally described. 

What is claimed as new and desired to be secured by 
Letters Patent of the United States is: 

1. A method of performing aeroelastic analyses, compris- 
ing the steps of: 

providing a mathematical model for a plurality of struc- 
tural modes of a flexible structure; 
providing a nonlinear aerodynamic model for generating a 
plurality of unsteady aerodynamic responses based on 
said plurality of structural modes for conditions defining 
an aerodynamic condition of said flexible structure; 
performing a single execution of said nonlinear aerody- 
namic model for said plurality of structural modes using 
a family of orthogonal functions as inputs thereto; 
converting said plurality of unsteady aerodynamic 
responses generated during said step of performing to a 
linear state- space model; and 
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generating static and dynamic aeroelastic solutions using 
computational interaction between said mathematical 
model and said linear state-space model for a plurality of 
periodic points in time. 

5 2. A method of performing aeroelastic analyses according 

to claim 1, wherein said generation of static and dynamic 
aeroelastic solutions comprises rigid aerodynamic offsets. 

3. A method of performing aeroelastic analyses according 
to claim 1, further comprising the steps of: 

10 changing a period for said plurality of periodic points in 
time; and 

regenerating static and dynamic aeroelastic solutions using 
computational interaction between said mathematical 

15 model and said linear state-space model for said plural- 
ity of periodic points in time so-changed in terms of said 
period thereof. 

4. A method of performing aeroelastic analyses according 
to claim 1, wherein said family of orthogonal functions is 

20 selected from the group consisting of block pulse functions, 
Haar functions, and Walsh functions. 

5. A method of performing aeroelastic analyses according 
to claim 1, wherein said nonlinear aerodynamic model is a 
computational fluid dynamics code. 

25 6. A method of performing aeroelastic analyses according 

to claim 1, wherein said step of converting is performed using 
an Eigensystem Realization Algorithm. 

7. A method of performing aeroelastic analyses according 
to claim 1, wherein said computational interaction comprises 

30 a feedback process. 

8. A method of performing aeroelastic analyses, compris- 
ing the steps of: 

providing a mathematical model for a plurality of struc- 

35 tural modes of a flexible structure; 

providing a nonlinear aerodynamic model for generating a 
plurality of unsteady aerodynamic responses based on 
said plurality of structural modes for conditions defining 
an aerodynamic condition of said flexible structure; 

40 generating a linear state-space model using a single execu- 
tion of said nonlinear aerodynamic model for said plu- 
rality of structural modes, said single execution of said 
nonlinear aerodynamic model using a family of orthogo- 
nal functions as inputs thereto; and 

45 generating static and dynamic aeroelastic solutions using 
computational interaction between said mathematical 
model and said linear state-space model for a plurality of 
periodic points in time. 

9. A method of performing aeroelastic analyses according 

50 to claim 8, wherein said generation of static and dynamic 

aeroelastic solutions comprises rigid aerodynamic offsets. 

10. A method of performing aeroelastic analyses according 
to claim 8, further comprising the steps of: 

changing a period for said plurality of periodic points in 

55 time; and 

regenerating static and dynamic aeroelastic solutions using 
computational interaction between said mathematical 
model and said linear state-space model for said plural- 
ity of periodic points in time so-changed in terms of said 

60 period thereof. 

11 . A method of performing aeroelastic analyses according 
to claim 8, wherein said family of orthogonal functions is 
selected from the group consisting of block pulse functions, 
Haar functions, and Walsh functions. 

65 12. A method of performing aeroelastic analyses according 

to claim 8, wherein said nonlinear aerodynamic model is a 
computational fluid dynamics code. 
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13 . A method of performing aeroelastic analyses according 
to claim 8, wherein said computational interaction comprises 
a feedback process. 

14. A method of performing aeroelastic analyses, compris- 
ing the steps of: 

providing a mathematical model for a plurality of struc- 
tural modes of a flexible structure; 

providing a nonlinear aerodynamic model for generating a 
plurality of unsteady aerodynamic responses based on 
said plurality of structural modes for conditions defining 
an aerodynamic condition of said flexible structure; 

performing a single execution of said nonlinear aerody- 
namic model for said plurality of structural modes using 
a family of orthogonal functions as inputs thereto, 
wherein said family of orthogonal functions is selected 
from the group consisting of block pulse functions, Haar 
functions, and Walsh functions; 

converting said plurality of unsteady aerodynamic 
responses generated during said step of performing to a 
linear state- space model; 

generating static and dynamic aeroelastic solutions using 
computational feedback interaction between said math- 
ematical model and said linear state-space model for a 
plurality of periodic points in time; 

changing a period for said plurality of periodic points in 
time; and 

regenerating static and dynamic aeroelastic solutions using 
computational feedback interaction between said math- 
ematical model and said linear state-space model for 
said plurality of periodic points in time so-changed in 
terms of said period thereof. 

15. A method of performing aeroelastic analyses according 
to claim 14, wherein said generation and regeneration of 
static and dynamic aeroelastic solutions comprises rigid aero- 
dynamic offsets. 

16. A method of performing aeroelastic analyses according 
to claim 14, wherein said nonlinear aerodynamic model is a 
computational fluid dynamics code. 

17. A method of performing aeroelastic analyses according 
to claim 14, wherein said step of converting is performed 
using an Eigensystem Realization Algorithm. 

18. A method of performing aeroelastic analyses, compris- 
ing the steps of: 

providing a mathematical model for a plurality of struc- 
tural modes of a flexible structure; 

providing a nonlinear aerodynamic model for generating a 
plurality of unsteady aerodynamic responses based on 
said plurality of structural modes for conditions defining 
an aerodynamic condition of said flexible structure; 

generating a linear state-space model using a single execu- 
tion of said nonlinear aerodynamic model for said plu- 
rality of structural modes, said single execution of said 
nonlinear aerodynamic model using a family of orthogo- 
nal functions as inputs thereto, wherein said family of 
orthogonal functions is selected from the group consist- 
ing of block pulse functions, Haar functions, and Walsh 
functions; 

generating static and dynamic aeroelastic solutions using 
computational feedback interaction between said math- 
ematical model and said linear state-space model for a 
plurality of periodic points in time; 
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changing a period for said plurality of periodic points in 
time; and 

regenerating static and dynamic aeroelastic solutions using 
computational feedback interaction between said math- 
5 ematical model and said linear state- space model for 
said plurality of periodic points in time so-changed in 
terms of said period thereof. 

19. A method of perfoming aeroelastic analyses according 
to claim 18, wherein said generation and regeneration of 

l o static and dynamic aeroelastic solutions comprises rigid aero - 
dynamic offsets. 

20. A method of performing aeroelastic analyses according 
to claim 18, wherein said nonlinear aerodynamic model is a 
computational fluid dynamics code. 

15 21 . A method of performing aeroelastic analyses, compris- 

ing the steps of: 

providing a state-space mathematical model for a plurality 
of structural modes of a flexible structure; 

providing a nonlinear aerodynamic model for generating a 
20 plurality of unsteady aerodynamic responses based on 
said plurality of structural modes for conditions defining 
an aerodynamic condition of said flexible structure; 

performing a single execution of said nonlinear aerody- 
namic model for said plurality of structural modes using 
25 a family of orthogonal functions as inputs thereto; 

converting said plurality of unsteady aerodynamic 
responses generated during said step of performing to a 
linear state-space model; and 

generating state-space static and dynamic aeroelastic solu- 
30 tions using computational interaction between said 
state-space mathematical model and said linear state- 
space model for a plurality of periodic points in time. 

22. A method of performing aeroelastic analyses according 
to claim 21, wherein said generation of state-space static and 

35 dynamic aeroelastic solutions comprises rigid aerodynamic 
offsets. 

23 . A method of performing aeroelastic analyses according 
to claim 21, further comprising the steps of: 

changing a period for said plurality of periodic points in 
40 time; and 

regenerating state-space static and dynamic aeroelastic 
solutions using computational interaction between said 
state-space mathematical model and said linear state- 
space model for said plurality of periodic points in time 
45 so-changed in terms of said period thereof. 

24. A method of performing aeroelastic analyses according 
to claim 21, wherein said family of orthogonal functions is 
selected from the group consisting of block pulse functions, 
Haar functions, and Walsh functions. 

50 25 . A method of performing aeroelastic analyses according 

to claim 21, wherein said nonlinear aerodynamic model is a 
computational fluid dynamics code. 

26. A method of performing aeroelastic analyses according 
to claim 21, wherein said step of converting is performed 

55 using an Eigensystem Realization Algorithm. 

27. A method of performing aeroelastic analyses according 
to claim 21, wherein said computational interaction com- 
prises a feedback process. 

60 



